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Inner vacuum gap formation and drifting subpulses in 
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Abstract. The problem of formation of the inner vaeuum 
gap in neutron stars with • B < is considered. It is 
argued by means of the condition Ti/Tg > 1, where Ti 
is the melting temperature of surface ions and Ts is 
the actual temperature of the polar cap surface, that the 
inner vacuum gap can form, provided that the actual sur- 
face magnetic field is extremaly strong {Bs> 10^"^ G) and 
curved [TZ < 10^ cm), irrespective of the value of dipo- 
lar component measured from the pulsar spin down rate. 
The calculations are carried out for pulsars with drifting 
subpulses and/or periodic intensity modulations, in which 
the existence of the quasi steady vacuum gap discharging 
via E X B drifting sparks is almost unavoidable. 



1. Introduction 

The subpulses in single pulses of a number of pulsars 
change phase systematically between adjacent pulses, 
forming apparent driftbands of the duration from sev- 
eral to few tenths of pulse periods. The subpulse intensity 
is also systematically modulated along driftbands, typi- 
cally increasing towards the pulse centre. In some pul- 
sars, which can be associated with the central cut of the 
line-of-sight trajectory, only periodic intensity modula- 
tion is observed, without a systematic change of subpulse 
phase. On the other hand, the clear subpulse driftbands 
are found in pulsars with grazing line-of-sight trajecto- 
ries. These characteristics strongly suggest an interpreta- 
tion of this phenomenon as a system of subpulse asso- 
ciated beams rotating slowly around the magnetic axis. 
Ruderman & Sutherland (1975; hereafter RS75) proposed 
a natural explanation of the origin of subpulse drift, that 
involved a number of isolated E x B drifting sparks dis- 
charging the quasi steady vacuum gap formed above the 
polar cap of the so-called anti-pulsar with Jl • B < 0, 
in which jeFe ions were strongly bound at the surface. 
Although the original idea of RS75 associating the rotat- 
ing sub-beams with the circulating sparks is still regarded 
as the best qualitative model of drifting subpulse phe- 
nomenon, their vacuum gap was later demonstrated to 
suffer from the so-called binding energy problem (for re- 



view see Abrahams & Shapiro 1991). In fact, the cohesive 



energies of"2gFe ions used by RS75 proved to be largely 
overestimated and the inner vacuum gap envisioned by 
RS75 was impossible to form. However, it is worth to em- 
phasize that RS75 considered the canonical surface dipo- 
lar magnetic fields with values determined from the pul- 
sar spindown rate, although they implicitly assumed small 
radii of curvature TZ ~ 10^ cm, inconsistent with purely 
dipolar field. Recently Gil & Mitra (2001; hereafter GMOl) 
revisited the binding energy problem an found that for- 
mation of vacuum gap is in principle possible, although 
it requires extremely strong non dipolar surface magnetic 
field Bs = bBd; where the coefficient 6 >> 1 in typical 
pulsar, Bd = 6.4 x 10^^{PP)° '^G = 2 x lO^^iPP^i^ f -^ G 
is the dipolar field at the pole, P is the pulsar period in 
seconds, P is the period derivative and P-15 = P/10~^^. 

In a superstrong surface magnetic field Bg > 0-lBq, 
where Bg = 4.414 x 10^"^ G, the asymptotic approximation 



of Erber 1966 used by RS75 in derivation of the height of 



quasi steady vacuum gap is no longer valid. In fact, in such 
strong field the high energy Ef = huj photons produce 
electron-positron pairs at or near the kinematic threshold 
huj = 2mc^/ sin6', where sin9 = h/TZ, h is the gap height, 
and TZ = T^-elO^ cm is the radius of curvature of surface 



magnetic field lines (e.g. Daugherty & Harding 1983), h is 
the Planck constant, c is the speed of light, m and e are the 
electron mass and charge, respectively. The vacuum gap 
formed under such conditions was called the Near Thresh- 
old Vacuum Gap (hereafter NTVG) by GMOl. They con- 
sidered two kinds of high energy seed photons dominat- 
ing the e^e+ pair production: the Curvature Radiation 
(CR) photons with energy hu = {3/2)n-f^c/TZ (RS75; 
Zhang fc Qiao 19*9^ ), and resonant Inverse Compton Scat- 
tering (ICS) photons with energy Huj = 2jheBs/mc (?), 
where 7 is a typical Lorentz factor of particles within 
the gap. The corresponding vacuum gap is called the 
Curvature Radiation dominated (CR-NTVG) and the In- 
verse Compton Scattering dominated (ICS-NTVG), re- 
spectively. GMOl estimated the characteristic heights of 
both CR-NTVG and ICS-NTVG. In this paper we will 
further refine these estimates by including the general rel- 
ativistic (GR) effects of inertial frame dragging (IFD) and 
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considering the heat flow conditions within the thin up- 
permost surface layer of the polar cap. Moreover, we will 
use broader range of cohesive energies of surface jgFe ions. 
The obtained VG models are applied to pulsars with drift- 
ing subpulses and/or periodic intensity modulations, in 
which the presence of E x B drifting spark discharges 
seems almost unavoidable (Deshpande & Rankin 1999; 
2001; Vivekanand & Joshi 1999 and Gil & Scndyk 2000). 

2. Near threshold vacuum gap formation 

If the cohesive energy of IgFe ions is large enough to pre- 
vent them from thermionic (this section) or field emission 
(Appendix A), a vacuum gap forms right above the polar 
cap. Jcssncr ct al. (2001) pointed out that the GR effect of 



IFD (?) should affect the RS75 type models with a vacuum 



gap above the polar cap. Although Jcssncr ct al. (2001) 



did not investigate the problem, they implicitly suggested 
that the electric fields distorted by GR effects make forma- 
tion of "starved" inner magnetospheric regions even more 
difficult than in the ffat space case. However, Zhang et al. 
(2000) demonstrated that although the GR-IFD effect is 
small, it nevertheless slightly helps formation of VG above 
the polar caps. In other words, the GR modified potential 
drop within the VG is slightly lower than in the flat space 
case. Below we conflrm this finding for NTVG conditions, 
with a very strong and complicated surface magnetic field. 

The gap electric field Ej| along Bg results from a de- 
viation of the local charge density p w from the co- 
rotational charge density pc = {(^/ajpcj, where poj — 
— r2 • Bs/(27rc) is the fiat space-time Goldreich- Julian 
(1969) charge density, ft ~ 27r/P, = 1 — Kg, Kg ^ 
{rg/R){I/MR^), a = (1 - rg/Rf/^ is the redshift factor, 
Tg is the gravitational radius, M is the neutron star mass 
and / is the neutron star moment of inertia. The potential 
V and electric fields E\\ within a gap are determined by 



GR-analog (Muslimov & Tsygan 1992) of the one dimen- 
sional Poisson's equation {1 / a)d'''V / dz^ = — 47r(p — pc) = 
2VlB sC, / {ca) , with a boundary condition E\\{z ^ h) = Q, 
where h is the height of infinitesimal gap. The solution of 
Poisson's equation gives 



and 



(1) 



(2) 



In further calculations we will adopt a typical value of the 
correction factor C ~ 0.85 (corresponding to M = IMq, 
R = 10^ cm and / = 10"'^ g cm^), although its value can 
be as low as 0.73 ( Zhang et al. 2000| ). Thus, the potential 



drop within the actual gap can be 15 to 27 percent lower 
than in the conventional RS75 model. 

The polar cap surface of the anti-pulsar (Jl • B < 0) is 
heated by a back-flow of relativistic electrons (accelerated 
in the parallel electric field i?||) to the temperature Tg = 



k^^'^{eAVN/anrpy^^, where AV is described by equa- 
tion (2), N = ■nr^Bs/eP is the Goldreich- Julian kinematic 
flux and the heat flow coefficient 0.2< fc < 1 is described 
in the Appendix B. The thermal condition for the vacuum 
gap formation can be written in the form Ti/Tg > 1, where 
Ts is the actual surface temperature described above, and 
Tj = A£c/30k is the iron critical (melting) temperature 
above which 26 Fe ions are not bound on the surface (?), 
where Aec is the cohesive energy of condensed 26Fe matter 
in the neutron star surface and k= 1.38 x 10~^^JK~^ is 
the Boltzman constant. The properties of condensed mat- 
ter in very strong magnetic flclds characteristic for the 
neutron star surface have been investigated by many au- 
thors using different e xamination methods (for review see 
Usov fc Melrose 1995 ). There exists a lot of discrepancy in 
determination of the cohesive energy Asc and in this paper 
we refer to the two most useful papers, representing the 
limiting extreme cases. Abrahams fc Shapiro (1991; AS91 
henceforth) estimated Asc = 0.91 keV, 2.9 keV and 4.9 
keV for Bs = lO^^ G, 5 x lO^^ G and 10^^ G, respectively. 



These values were approximated by Usov fc Melrose 1995 
in the form Ae^ ^ (0.9keV)(B,,/10i2 Gf-'^^, which leads 
to melting temperatures 

Ti = (3.5 X 10^K)(S,/1012G)°-^3 = 
(6 X 105)6°-^3(pp_^^)0.36 ^ 



On the other hand, Jones (1986; J86 henceforth) ob- 
tained much lower cohesive energies Aec = 0.29 keV, 
0.60 keV and 0.92 keV for = 2 x 10^^ g, 5 x lO^^ G 
and 10^^ G, respectively. This can be approximated by 
Aec ^ (0.18 keV)(_Bs/10i2Q)0.7 g^^j converted to melting 
temperatures 



l2r'^o.7 



T, = (0.7 X 10^K)(Bs/10^^G) 
(1.2 X 105)6" '^(P- P-if 



\0.36 



(4) 



Below we consider the condition Ti/Tg > 1, using both 
expressions for melting temperatures described by equa- 
tions (3) and (4), for CR- and ICS-dominated NTVG mod- 
els, separately. 

2.1. CR-NTVG 

In this case the gap height h = hcR is determined by the 
condition that h = Iph, where Iph ~ sinO TZ = {B±/ Bs)TZ 
is the mean free path for pair production by photon prop- 
agating at an angle 6 to the local surface magnetic field 
(RS75). The CR-NTVG model is described by the follow- 
ing parameters: the height of quasi steady gap 



T-CR 



(3 X lo3)c-3/^7^^/'6-3/^p3/l^p^'/'' 



(5) 



(notice typographical errors in eq. [6] of GMOl, which are 
corrected here), the gap potential drop 

AV = (1.2 X 10i2)Ci/%/'ri/^p-i/i4pi/i4 
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and the surface temperature 

Ts = (3.4 X 106)Cl/28fcl/47^1/7^2/7p-l/7pl/7 ^ 

The thermal condition Ti/Tg > 1 for the formation of 
CR-NTVG leads to a family of critical lines on the P — P 
diagram (see Fig. 1 in GMOl) 

> ^'c°-i^fci-i^<^^fe-'p-' ■^ (8) 

where A = (2.7 x IQ^f^ = 52 for AS91 case (eq. [3]) 
and A = (3.96 x lO^)!/^ = 1990 for J86 case (eq. [4]). 
Alternatively, one can find a minimum required surface 
magnetic field Bg = bBd expressed by the coefficient b in 
the form 

= ^c°■°«fc°■"7^r'^-^■"^^'--°5'■ (9) 



2.2. ICS-NTVG 

In this case the gap height h ~ hjcs is determined by 
the condition h = Iph ~ le, where le is the mean free 
path of the electron to emit a photon with energy hu) = 



2jheBs/mc (GMOl, [Zhang ct al. 200q ). The ICS-NTVG 
model is described by the following parameters: the height 
of quasi steady gap 



hics = (5 x lQ'')C'k-^ '"n''e''b-'p-^-'''PI^r!' cm. 



the gap potential drop 
AV = (5.2 x 10l2)cO■^2^-o■l47^l•"5-l 

and the surface temperature 

6\a0.18 7„0.25x>0.28 



(10) 



= (4 X 10'')C'-'''fc 



(11) 



(12) 



The thermal condition Ti/Ts > 1 for the formation of 
ICS-NTVG leads to a family of critical lines on the P — P 
diagram (see Fig. 1 in GMOl) 



P-15 > B'C 



2^0.5^0.77^0.8^-.p- 



-2.2 



(13) 



where 6 = (2 x 10^)1/2 ^ 14 for AS91 case (eq. [3]) and 
B = (1.69 X 10**)i/2 ^ 130 for J86 case (eq. [4]). Alterna- 
tively, one can find a minimum required surface magnetic 
field Bs = bBd expressed by the coefficient 6 in the form 



0.5 



(14) 



2.3. NTVG development in pulsars with drifting 
subpulses 

GMOl examined the P — P diagram (their Fig. 1) with 
538 pulsars from the Pulsar Catalog (Taylor et al. 1993) 



with respect to the possibility of VG formation by means 
of the condition Ti/Tg, with Ti determined according to 
AS91 (eq. [3]). They concluded that formation of VG is 
in principle possible, although it requires very strong and 



curved surface magnetic field Bg = b-Bd> 10^^ G, irrespec- 
tive of the value of Bd. Here we reexamine this problem 
by means of equations (9) and (14), with melting temper- 
atures corresponding to both the AS91 case {A = 52 and 
B = 14) and J86 case {A = 1990 and B = 130), for CR 
and ICS seed photons, respectively. Therefore, we cover 
practically almost the whole range of melting tempera- 
tures between the two limiting cases, which differ by a 
factor of five between each other. We adopt the GR-IFD 
correction factor C = 0.85 (?), the normalized radius of 
curvature of surface magnetic field lines 0.01 < TZq < 1.0 
(GMOl and references therein) and the heat flow coeffi- 
cient 0.2 < fc < 1.0 (Appendix B). 

The results of calculations of NTVG models for 42 
pulsars with drifting subpulses and/or periodic intensity 
modulations (after Rankin 1986) are presented in Fig. 1. 
We calculated the ratio Bs/Bg = 0.0453 6(PP_15)°■^ 
where Bq = 4.414 x 10^'^ G and the coefficient b is de- 
termined by equation 9 and equation 14 for CR and ICS 
seed photons, respectively. Noting that the functional de- 
pendence on P and P-15 in both equations is almost iden- 
tical, we sorted the pulsar names (shown on the horizontal 
axes) according to the increasing value of the ratio Bg/Bg 
(shown on the vertical axes). Calculations were carried 
out for Bg/Bq > 0.1, since below this value the NT con- 
ditions are no longer valid. On the other hand, the values 
of Bg/Bq are limited from above by the photon splitting 
threshold, which is roughly about 10"^^ G (see also astro- 
ph/0102097, ?). Therefore, the physically plausible NTVG 
models lie within the grey areas, which correspond to sur- 
face magnetic fields 4.4 x 10^2 G < Bg < 10" G. The 
left hand side of Fig. 1 corresponds to AS91 case {A = 52 
and S = 14 in cqs. [9] and [14], respectively) and the right 
hand side of Fig. 1 corresponds to J86 case {A — 1930 
and B = 130 in eqs. [9] and [14], respectively). Four pan- 
els in each side of the figure correspond to different val- 
ues of the radius of curvature TZq = 1.0, 0.1. 0.05 and 
0.01 from top to bottom, respectively (indicated in left 
upper corner of each panel). Two sets of curved lines in 
each panel correspond to the CR (thin upper lines) and 
ICS (thick lower lines) seed photons, respectively. Three 
different lines within each set correspond to different val- 
ues of the heat flow coefficient k ~ 1.0 (dotted), k = 0.6 
(dashed) and k ~ 0.2 (long dashed). 

A visual inspection of the model curves within grey 
areas in Fig. 1 shows that in AS91 case the ICS-NTVG is 
favored for larger radii of curvature TZq > 0.05, while the 
CR-NTVG requires lower values of TZe < 0.1. In J86 case, 
only the ICS-NTVG corresponding to TZq < 0.1 can de- 
velop. If the actual cohesive energies correspond to some 
intermediate case between the AS91 and J86 cases, they 
will also be associated with the ICS-NTVG. Therefore, 
wc can generally conclude that the ICS-NTVG model 
is apparently favored in pulsars with drifting subpulses. 
The CR-NTVG is also possible although it requires an 
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extremely strong [Bs/Bq ~ 1) and/or extremely curved 
(7^6 ^ 0.01) surface magnetic field at the polar cap. 

3. Conclusions 

There is a growing evidence that the radio emission of pul- 
sars with systematically drifting subpulses (grazing cuts of 
the line-of-sight) or periodic intensity modulations (cen- 
tral cuts of the line-of-sight) is based on the inner vacuum 
gap developed just above the polar cap (?). To overcome 
the binding energy problem (?) put forward an attractive 
but exotic conjecture that pulsars showing the drifting 
subpulses represent bare polar cap strange stars (BPCSS) 
rather than neutron stars. However, as demonstrated in 
this paper, invoking the BPCSS conjecture is not neces- 
sary to explain drifting subpulse phenomenon. The quasi 
steady vacuum gap, with either curvature radiation or in- 
verse Compton scattering seed photons, can easily form 
in pulsars with Jl • B < 0, provided that the actual sur- 
face magnetic field at the polar cap is extremely strong 
Bs ^ 10^'^ G and curved TL < 10^ cm, irrespective of 
the value of dipolar component measured from the pulsar 
spindown rate. 

The pulsars with drifting subpulses and/or periodic 
intensity modulation do not seem to occupy any partic- 
ular region of the P — P diagram (see Fig. 1 in GMOl). 
Rather, they are spread uniformly all over the P—P space, 
at least for typical pulsars. Therefore, it seems tempt- 
ing to propose that radio emission of all pulsars is driven 
by vacuum gap activities. An attractive property of such 
proposition is that the nonstationary sparking discharges 
induce the two-steam instabilities that develop at rela- 
tively low altitudes (?), where the pulsars radio emission 
is expected to originate (?). It is generally believed that 
the high frequency plasma waves generated by the two- 
stream instabilities can be converted into coherent elec- 
tromagnetic radiation at pulsar radio wavelengths (e.g. 
Melikidze et al. 2000 ). In such a scenario, all radio pul- 
sars would require a strong, non-dipolar surface magnetic 
field at their polar caps (e.g. Gil et al. 2002 a, b). 
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dr B. Zhang for helpful discussions. We also thank E. Gil, U. 
Maciejewska, and M. Sendyk for technical help. 

Appendix A: Field emission 

The vacuum gap can form in pulsars with • B < 
if the actual surface temperature Tg is not high 
enough to liberate 26^6 ions from the polar cap sur- 
face by means of thermal emission. Now we will ex- 
amine the field (cold cathode) emission, which is pos- 
sibly important when thermionic emission is negligi- 
ble (e.g. [Zhang fc Harding 2000| ). The maximum elec- 
tric field (along Bg) at the NS surface i?|[(max) 

^6 



or taking into account that Q ~ 0.85 and b> 2(P-P_i5) 
E\\{max) < (5 X W)nl-^'^ P^-^^P'^f^ V/cm. The critical 
electric field needed to pull |gFe iron ions from the NS 
surface is E\\{crit) = (8 x 10i2)(Ae/26 keV)^/^ V/cm 
(?), where the cohesive (binding) energy of iron ions in 
strong surface magnetic field Bg ^ 10^^ G is Ae = 



4.85 keV ( [Abrahams fc Shapiro 199l| ). Thus, E\\{crit) = 
6.4 x 10^^ V/cm >> ~ 10''' V/cm and no field emission 
occurs. It is possible that the cohesive energy is largely 
overestimated and Aec can be much smaller than about 
4 keV even at Bs > 10^"^ G. We can thus ask about the 
minimum value of Asc at which the field emission is still 
negligible. By direct comparison of (max) and (crit) 
we obtain Ae^ > 40a;"-S'^ eV, where x = 7^g■5^pO■6■*p0.5. 
is of the order of unity. Therefore, contrary to the conclu- 
sion of lessner et al. (2001), no field emission is expected 
under any circumstances. 

Appendix B: Heat flow conditions at the polar 
cap surface 

Let us consider the heat flow conditions within the up- 
permost surface layer of the pulsar polar cap above which 
NTVG can develope. T he basic heat flow equation is (e.g. 
Eichlcr fc Cheng 1989[) 



C— ~ — fn — 

dt dx \ " dx 



(Bl) 



where C is the heat capacity (per unit volume) and 
K|| >> is the thermal conductivity along (||) sur- 

face magnetic field lines (which are assumed to be perpen- 
dicular (_L) to the polar cap surface). The heating of the 
surface layer of thickness Ax is due to sparking avalanche 
with a characteristic development time scale At. We can 
write approximately dT/dt « T/At, dT/dx « T/Ax, 
d/dx{dT/dx) « T/Ax^ and thus C/At « k\\Ax^. There- 
fore, the crust thickness Ax heated during At is approxi- 

1/2 

mately Ax « (K||Ai/C) (within an uncertainty factor 
5x of 2 or so), where the time scale corresponds to the 
spark development time At « 10/xs (RS75). The energy 
balance equation is Qheat = Qrad + Qcond, where Qheat = 
encjAVmax (e.g. RS75), Qrad = crT^ (Ts is the actual sur- 
face temperature and cr = 5.67 x 10~^ erg cm~'^K~'^s~^), 
and Qcond = —n\\dT/dx ~ —k\]Ts/ Ax. We can now define 
the heat flow coefficient k = Qrad/ Qheat < 1, which de- 
scribes deviations from the black-body conditions on the 
surface of the sparking polar cap. In other words, the value 
of k describes the amount of heat conducted beneath the 
polar cap which cannot be transferred back to the surface 
from the penetration depth Aa; during the time-scale At. 
The coefficient k can be written in the form 



1 



1 



(:{4Tr/cP)Bghics = (1.25 x 109)C°■866-l7^g■"pO.l4 y/cm. 



l-|-K||/(crT3Aa 



1-f (K||C)l/2/(aT3Atl/2)' 

(B2) 
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Surface temperature 

Fig. Bl. Heat flow coefficient k versus surface tempera- 
ture Tq (see text for explanations). 



whose value can be estimated once the parameters C, K|| 
and Ts as well as Aa; or At are known. 

The matter density at the neutron star surface 
composed mainly of 2|Fe ions (e.g. |Usov fc Melrose 1995 ) 
i s p{Bs) - 4 X 10^ {b7/T(P G)^/5g cm-3 
( [Flowers et al. 19771 ). Thus for ~ (1 ^ 3) x 10" G 
we have p ~ (0.6 2.4) x lO^g cm^^. The thermal 

2.2 X 

' 5 



energy density of |j?Fe is Upe — 2.2 x lO^^p^Ts erg cm^^ 
( Eichlcr fc Cheng 1989 ), where p5 = p/(10^g cm~^) and 

dUpe/dt = 
10^^P5 erg cm~'^K~^. For 



Ts = r,/(10^ K). Thus, the heat capacity C 
IQ-^dEpe/dTs ^ (1 ^ 5) X iriii- 
ps w 1 we have C w 2 X lO^^erg cm-^K'^ The 
longitudinal thermal conductivity can be estimated 
as K\\ (2-^4) X 10^^ erg s^^cm^^K^^ (see Fig. 5 
in [Potekhin 1999| ). Thus K||/C « 1.5 x 10^ cm^ and 
the penetration depth Ax « lOAt^/^s^^/^ cm. Since 
At ~ 10/i = 10^^ s (spark characteristic time scale) then 
Aa; « 0.03 cm. More generally, the penetration depth 
can be written as 



Ax ^ 0.03 Sx 



«13 



1/2 



(B3) 



where K13 = k^/IO^^ w 2 ^ 4 and Cn = « 1 ^ 5 

and the uncertainty factor Sx « 0.5 -f- 2. Now the heat flow 
coefficient can be written as 



from independent considerations that in PSR B0943-I-10 
the heat flow coefficient k < 0.8, in consistency with the 
results obtained in this paper. 
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where Tg = T,/10^ and A = (KiaCnj^/Vfe « 0.7 ~ 6.3. 
Figure 2 shows variations of the heat flow coefficient k ver- 
sus the surface temperature Tq (in 10^ K) for three values 
of A = 1 (upper curve), A = 3 (middle curve) and A = 6 
(lower curve). As one can see, for realistic surface temper- 
atures of few times 10^ K, the value s of the heat flow coef- 
flcient k are in the range 0.2 0.9. Gil et al. 2002a foimd 



